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Massive stars produce copious quantities of ultraviolet radiation beyond the Lyman limit, pho- 
toionizing the interstellar medium (ISM) and producing H II regions. As strong sources of 
recombination- and forbidden-line emission, infrared continuum, and thermal (free-free) radio 
continuum, H II regions serve as readily-observable beacons of massive star formation in the 
Milky Way and external galaxies. Along with supemovae, H II regions are dominant sources of 
^ feedback in star-forming galaxies, injecting radiative and mechanical luminosity into the ISM. 

H II regions may prove more important than supernovae as triggers of star formation through lo- 
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calized compression of cold cloud cores. In this review, I give a broad overview of the structure 
and time-evolution of H II regions, emphasizing complications to the theoretical picture revealed 
by multiwavelength observations. I discuss a recent controversy surrounding the dominant feed- 
back mechanism in 30 Doradus, the most luminous H II region in the Local Group. I summarize 
the first results from the Milky Way Project (MWP), which has produced a new catalog of several 
thousand candidate Galactic H II regions by enlisting >35,000 "citizen scientists" to search Spitzer 
Space Telescope survey images for bubble-shaped structures. The MWP and similar large cata- 



o 



f--->) logs enable empirical studies of Galactic H II region evolution across the full range of luminosities 

and statistical studies of triggered star formation. 
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1. Introduction: The Interstellar Medium Out of Balance 

Most students of astronomy will encounter a course on the physics of the interstellar medium 
(ISM; the gas and dust occupying the space between the stars) at some point in their undergraduate 
or post-graduate studies. Hence the concept that the ISM exists in several distinct "phases," in 
rough pressure equilibrium with each other, may be familiar to most readers. Physical quantities 
characterizing the major phases of ISM gas in the MiUcy Way are summarized in Table 1, compiled 
from the textbooks by Tielens [1] and Draine [2]. The coexistence of gas at such dramatically 
different densities, temperatures, and ionization fractions in the same galaxy is explained by the 
different filling factors of each phase. While the precise values for the filling factors remain poorly- 
measured and controversial, the basic picture that the colder, denser phases of the ISM exists as 
smaller clouds within the more diffuse, warm/hot phases, is well established. H II regions, localized 
regions of photo-ionized gas produced by hot, massive, OB-type stars, occupy a negligible fraction 
of the ISM volume, hence perhaps H II regions ought not to be regarded as a proper ISM phase 
at all. However, massive stars form in the densest regions of cold, molecular clouds, and as their 
far-ultraviolet (UV) radiation first photo-dissociates molecules and then photo-ionizes atoms, the 
multi-phase physics of the ISM can be studied within a single, small volume. 



Phase 


Density 


T 


Total Mass 


Scaleheight 


Filling 




(cm^^) 


(K) 


(10" Mq) 


(pc) 


factor 


Hot ionized medium 


~0.004 


~10^ 




3000 


~0.5 


Warm neutral medium 


0.5-0.6 


8000 


2.8 


~300 


~0.4 


Warm ionized medium 


0.1-0.3 


~5000 


1.0 


900 


~0.1 


Cold neutral medium 


30-50 


80-100 


2.2 


100 


~0.01 


Molecular Clouds 


10^-10^ 


10-50 


1.3 


75 


-10-4 


H II Regions 


1-104 


104 


0.05 


70 





Table 1: Phases of ISM Gas in the Milky Way 



In one of the immortalized insights of early modern astrophysics, Stromgren [3] realized that 
the mean free path in neutral hydrogen of UV photons beyond the Lyman limit (A < 912 A) is negli- 
gibly small compared to the size of the ionized hydrogen region produced by a hot star. H II regions 
therefore have sharp boundaries, or I-fronts, where the recombination rate balances the ionization 
rate. For a single star in an ambient medium of constant density, this boundary is defined as the 
Stromgren radius, 

f 3(2o \ 

where Qq is the ionizing photon rate (which depends on the star or stars responsible for the H II re- 
gion), nu is the hydrogen gas density, and qb is the Case B recombination coefficient [4]. Generally 
speaking, main-sequence or giant stars earlier than B3 emit sufficient to produce observable 
Galactic H ll regions. 

Because the ionizing stars provide an internal source of radiative and mechanical luminosity, 
H II regions rapidly become overpressured compared to the ambient ISM and expand. If the dom- 
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inant source of pressure is coUisional heating of gas by free electrons, then the time-evolution of 
the expansion follows the simple analytic relation from Lyman Spitzer's classic ISM text [4]: 
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where Cs2 ~ 10 km s^^ is the sound speed in the ionized gas. The expansion velocity obtained 
by differentiating this relation can exceed the (significantly lower) sound speed in the ambient 
medium, hence expanding I-fronts often become shock waves. If the ambient medium is molecular, 
UV photons emerging from the H II region destroy the molecules, creating a photodissociation 
region (PDR) around the H II region. 

Early O stars and OB giants drive powerful winds that fundamentally alter the structure of 
more luminous H II regions. Castor et al. [5] and Weaver et al. [6] provided analytical models for 
wind-blown bubbles produced by isolated, hot stars, and numerous authors have subsequently re- 
fined these models using a variety of semi-analytical and numerical techniques [7-11]. The basic, 
"onion- layer" structure of a wind-blown bubble is illustrated in Figure 1 . The highly supersonic 
(1000-2000 km s^') stellar wind flows freely outward for a short distance Rw from the star before 
it is shocked, producing a bubble of very hot, ionized gas. The "classical" H II region collapses into 
a photoionized shell of gas, (imperfectly) separated from the hot gas zone by a contact discontinuity 
at Re- The I-front at Rjp still represents the outer boundary of the H II region. 
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Outer, I-front shock (R;f) 



SHOCKED 
INTERSTELLAR GAS 
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Hot! >I07 K; X-rays 
Wind shock (Rw) 

Not to scald 



Figure 1: Anatomy of a wind-blown bubble, adapted from Weaver et al. [6] 



Many complications separate the ideal models of Stromgren spheres and wind-blown bubbles 
from reality. The ISM is clumpy, hence the ambient medium surrounding an H II region is never 
uniform. Massive stars tend to form in clusters, hence multiple stai^s often contribute to the ioniza- 
tion of a single H II region. The ambient ISM is generally in motion with respect to the ionizing 
star(s), and stellar winds need not be spherically symmetric. Dust mixed with gas in H II regions 
enables dramatic, radiative cooling [12], and magnetic fields threading through the clouds con- 
tribute anisotropic pressure support [13]. Turbulence provides additional pressure and facilitates 
mixing at the interfaces between gas layers. Early models neglected completely the contribution of 
radiation pressure [14, 15]. 
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Bubble N49 




Figure 2: Multiwavelength images of the wind-blown bubble N49 [17] and the giant H II region M17 [18], 
displayed at approximately the same physical scale. N49 image: red = SpitzerfMWSi 24 /im, green/blue = 
Spitzer/JRAC 8.0/4.5 nm. Mil image: red = MSX 21.3 /zm, green = IRAC 5.8 jUm, blue = Chandra soft 
(0.5-2 keV), diffuse X-rays. In both images, diamonds denote known O and early B stars and contours show 
20 cm thermal radio continuum. 



2. IMultiwavelength Observations of H ii Regions 

In spite of the messy complexity governing the structure of real H II regions, the basic struc- 
tures predicted by the wind-blown bubble models are identifiable in modern, multiwavelength 
images. In Figure 2, a prototypical wind-blown bubble ionized by a single 06.5 V star (Qq = 
8.5 X 10^^ s^^) [16, 17] is compared to the giant H II region M17, ionized by a dozen O stars, 
including several 04 V stars (Qq = 3x 10^° s"') [18]. In N49, the I-front at Rjp is defined by the 
sharp inner rim of 8.0 /im (green) polycyclic aromatic hydrocarbon (PAH) emission from the PDR, 
neatly encapsulating the photoionized gas shell (contours). Dust mixed with the photoionized gas 
and heated by radiation from the central star forms a torus of 24 jUm emission. Everett & Church- 
well [12] found that the lifetime of dust grains within the harsh environment of this H II region is 
extremely short, and suggested that dust must be continuously replenished from evaporating dense 
clumps to produce the observed 24 jxra emission. Draine [15] demonstrated that radiation pressure 
can produce evacuated cavities in the centers of H II regions, but noted that the central hole in N49 
is too large to be explained by radiation pressure alone and suggested that the stellar wind also 
contributes. An interface analogous to the contact discontinuity Rc in Figure 1 is likely located 
within the radio shell/24 jUm torus in N49, hence Rc/Rif < 1/2. 

Unlike N49, M17 is far from round, yet similar moiphological features can be discerned, 
with one important addition (Figure 2). A spectacular plume of hot. X-ray-emitting plasma (blue) 
occupies the central cavity of Ml 7 [19]. This X-ray emission provides direct evidence for stellar 
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wind shocks. At an absorption-corrected X-ray luminosity Lx = 1 x lO^*^ erg s^^ this plasma 
is unusually bright in comparison to other H II regions, but it is fainter than the predictions of 
wind-blown-bubble theory by more than an order of magnitude [20]. This discrepancy may be 
explained by collisional interactions with dust grains providing a mechanism for cooling the hot 
plasma, and/or depressurization of the wind-blown bubble where the plasma is not completely 
confined by the nebula. Either interpretation implies that the contact discontinuity in Figure 1 
does not effectively separate the photoionzed gas/dusty shell from the hot gas bubble in M17. 
Again assuming Rc falls at the inner edge of the photoionized shell and heated dust emission, 
Rc/^iF ~ 1/2 in M17. The inner cavity in M17 is clearly larger (both in absolute volume and as a 
fraction of the H II region volume) in comparison to that of N49. The photoionized shell in Ml 7 is 
supported by a combination of radiation pressure and hot gas pressure [13]. 

Although N49 and M17 are representative of a range of Galactic H II regions where stellar 
winds play an important role, on the Galactic scale feedback is dominated by the most luminous, 
starburst regions. In Figure 3, a Galactic starburst region, W43, is compared to 30 Doradus in the 
Large Magellanic Cloud (LMC), the most luminous H II region in the Local Group {Qq = 4.2 x 10^^ 
s^'; [21]). Mid-infrared (IR) images of the starburst H II regions reveal layered bubble morpholo- 
gies that are remarkably similar between the two regions. The large bubble lobes shown in each 
panel of Figure 3 are part of larger H II region complexes, with the ionizing clusters partially (in 30 
Dor) or completely (W43) obscured by dense, foreground filaments of bright, mid-IR emission. In 
the 30 Dor image, the PDRs appear pink and the photo-ionized shells green, while the W43 image 
matches the color-code of Figure 2, in which the PDRs appear yellow-green and dust mixed within 
photoionized shell appears red. Here it is most appropriate to describe the photoionized gas struc- 




Galactic Starburst 
Region W43 



1 



(Not to scale) 

Figure 3: Spitzer mid-IR images of two starburst H II regions. Left: Image of 30 Dor in the Large Magellanic 
Cloud, red = 8.0 jim, green = 4.5 ^m, blue = 3.6 jim. Right: Image of W43, red = 24 /im, green = 8.0 lJ.m, 
blue = 4.5 lim. 
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tures as shells, occupying thin layers just interior to the PDRs, with Rc ~ Rip. Indeed, the bubbles 
in 30 Dor are known to be filled with hot. X-ray-emitting plasma [22], and similar plasma would 
likely be found in W43 if comparable Chandra X-ray Observatory observations were obtained. 

Comparing Figures 2 and 3, a trend becomes apparent: as the luminosity of the ionizing cluster 
increases, so does Rc/Rif, the size of the central cavity relative to the overall size of the H II region. 
We may extend this trend down to low-luminosity H II regions ionized by late O or early B stars, 
for which the central cavity disappears entirely (Rc = 0). It will be useful to bear this trend in mind 
when considering the controversial issue of precisely which feedback mechanism, hot gas pressure 
or direct radiation pressure, dominates H II region structure at the high-luminosity extreme. 

3. A Rumble in the Tarantula: What is the Dominant Feedback Mechanism in 30 
Doradus? 

Thanks to its status as the most luminous H II region in the Local Group and its location in the 
low-metallicity environment of the LMC, 30 Dor (popularly known as the Tarantula Nebula) has 
long received intense observational scrutiny. At J = 50 kpc, 30 Dor is the best nearby laboratory 
for studying the physical conditions that prevailed in the unresolvable, high-redshift star-forming 
regions that dominated the major cosmological epoch of galaxy-building [21]. Recently, Lopez 
et al. [23, hereafter Lll] and Pellegrini et al. [21, hereafter Pll] canied out independent, pai^allel 
studies of the feedback processes shaping 30 Dor. Using fundamentally different approaches to 
interpreting multiwavelength datasets, these authors reached diametrically opposed conclusions; 
Lll reported that direct stellar radiation pressure dominates the interior of the H II region, while 
PI 1 argued that the pressure of the hot. X-ray-emitting plasma shapes the large-scale structure and 
dynamics. This disagreement is rooted in the different definitions of radiation pressure and the 
different assumed nebular geometries used in the two studies. 

Lll used the simplest definition of direct radiation pressure, 

Pd. = l/^ (3.1) 

(their equation 1), where Lboi is the bolometric luminosity of each star and r is the distance traveled 
by the starUght to reach a given point in the nebula, deprojected assuming a spherical geometry. Pjir 
declines sharply with distance from the central star cluster, R136 (note that this expression diverges 
for r = 0). 

By contrast, Pll constructed a non-symmetric, cavity model (based on the central region of 
30 Dor shown in Figure 3) for the nebular geometry and used photoionization models to calculate 
the density of H atoms nn and hence the ionization parameter U at each position, 

ATlr'^criH 

The divergent behavior of this expression is avoided by implementing the cavity model, in which 
the ionized gas is confined to shell structures near the I-fronts, and U is not calculated for the cavity 
interiors, where r^nn — )• 0. Pll then approximated the pressure exerted on the observed ionized 
gas by starlight in terms of the ionization parameter as 

P.i.r. = UnH{hv)^, (3.3) 
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where Lboi and Lq are the total bolometric and ionizing photon luminosity for all stars (assumed to 
be centered at R136) and (/zv) ~ 20 eV is average energy per ionizing photon (their equation 8). 

Ahhough Equations 3.3 and 3.2 can be combined and trivially reduced to Equation 3.1, doing 
so hides the ambiguous role of radiation pressure in regions where he vanishes. Pstars as definited 
by Pll represents the momentum imparted to the observed nebular gas. Lll acknowledge this 
alternative definition of radiation pressure, but claim that "it is necessary to characterize fdir as the 
energy density of the radiation field, since that definition reflects the total energy and momentum 
budget available to drive motion." This definition implies that the luminosity emitted by the OB 
stars could impart momentum with 100% efficiency everywhere in the nebula at once, an ideal case 
that could never occur in a real H II region. The justification brings to mind the old philosophical 
thought experiment about whether a tree falling in a forest makes a sound if there is no one around 
to hear it. If one dropped a cloud of dense, neutral gas close to the R136 star cluster, it would 
experience an enormous radiation pressure. But there are no dense clouds of neutral gas within 
the radiation-pressure-dominated region identified by LI 1. Instead, the interior structure of 30 Dor 
consists of evacuated cavities filled with hot, highly ionized plasma [22], and pressure from either 
the hot. X-ray plasma or the warm, photoionized gas can exceed the radiation pressure on the cavity 
walls (Lll, PI 1). Radiation pressure could have dominated the expansion of these large cavities in 
the past, when they were smaller (PI 1). It is difficult to answer this question definitively because 
the hot gas pressure in 30 Dor remains uncertain. 

Both Lll and Pll reported values for the hot gas pressure in 30 Dor. The pressure of the 
X-ray-emitting plasma can be calculated as 

Px = l.9nxkTx, (3.4) 

where nx and Tx are the density and temperature of the X-ray-emitting plasma. Pll used the 
results of spectral fits by Townsley et al. [22] while Lll performed their own spectral fitting to 
the same archival data plus a newer, 90-ks Chandra observation of 30 Dor [24]. It is difficult to 
derive nx accurately from the emission measure returned by spectral fitting because it depends 
strongly on the assumed geometry of the X-ray bubbles. Lll treated 30 Dor as a "beach ball" 
with a global, spherical geometry, which effectively minimizes nx- PI 1 treated 30 Dor as a "bunch 
of grapes," assuming a spherical geometry for each smaller, diffuse X-ray region identified by 
Townsley et al. [22], which yields higher nx- Spectral modeling of diffuse emission structures in 
regions like 30 Dor is a tremendously complicated task [20], as any given sightline will contain 
plasma from multiple origins, including both stellar wind shocks and supernova, at a variety of 
temperatures, densities, and compositions. To illustrate the pitfalls of over-interpreting these data, 
I note that different approaches to the spectral fitting of the the global, diffuse X-ray emission 
of 30 Dor yield significantly different results. Lll fit a single-temperature plasma model and 
reported kTx = 0.64^q q2 keV and absorption-corrected Lx = 4.5 x 10^^ erg s^'. In contrast, the 
most recent spectral fits by L. K. Townsley (private communication) employ 3 plasma components 
(plus numerous gaussian profiles to fit unidentified emission lines) ranging from 0.3-0.8 keV and 
yield Lx = 1.2-1.9 x 10^' erg s^^. It is particulai^ly difficult to discern whether regions of 30 
Dor that appear X-ray dark represent the boundaries of confined, hot plasma bubbles, or whether 
the plasma extends behind foreground material that absorbs the soft X-rays. The existing 114 ks 
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combined Chandra integration represents a very shallow observation at the distance of 30 Dor 
when compared to observations of Galactic H II regions. 

The story of 30 Dor has a moral; given the complexity and ambiguity involved in interpreting 
the multiwavelength data on this resolved starburst region, investigators wishing to extend these re- 
sults to draw conclusions about feedback mechanisms shaping unresolved regions at cosmological 
distances do so at their own risk. 

4. The MUky Way Project H ii Region Catalog 

Strong empirical constraints on the time-evolution of feedback-driven H II regions require the 
comparative study of large observational samples of H II regions. The wealth of new IR imaging 
data provided by the Spitzer Galactic Legacy Infrared Mid-Plane Extraordinaire (GLIMPSE) [25] 
and subsequent high-resolution surveys of the Milky Way allow us to penetrate the obscuring veil 
of dust in the Galactic plane, revealing the structure of H II regions and PDRs in unprecedented 
detail. Using the GLIMPSE images, Churchwell et al. [16, 26] cataloged neai^ly 600 IR bubbles, 
ring and arc-shaped structures apparent in 8 /im PAH emission. The majority of these structures 
ai^e PDRs surrounding H II regions, from energetic, wind-blown bubbles to low-luminosity nebulae 
surrounding B-type stars. 

Over the past year, >35,000 internet users from around the world have been searching for 
bubbles in Spitzer survey images of the Galactic plane as part of the Milky Way Project (MWP; 
http://www.milkywayproject.org), a recent installment in the Zooniverse, the premier series of on- 
Une "citizen science" projects (http://www.zooniverze.org). Upon creating a Zooniverse account 
and logging into milkywayproject.org, MWP volunteers are presented with a random image and 
instructed to identify and fit structures within that image that resemble bubble rims (PDRs resem- 
bling the regions in Figure 2) with elliptical annuli (or mark the locations of bubbles that are barely 
resolved with boxes). By carefully combining the results from many individuals for each part of 
the sky, the MWP simultaneously leverages the superior pattern-recognition skills of the human 
eye-brain system and takes advantage of the "wisdom of crowds." The MWP First Data Release 
presents a catalog of 5,106 H II regions, representing an order of magnitude improvement in com- 
pleteness compared to previous catalogs [27]. 

The MWP bubbles catalog will facilitate the study of triggered star formation on the Galactic 
scale. The idea that feedback from expanding H II regions can exert external pressure on cold, 
molecular cloud cores, resulting in self-propagating, sequential massive star formation has been 
around for decades [28-31] but recently has seen a resurgence of observational and theoretical 
interest, motivated in large part by the identification of numerous instances of small bubbles, young 
stellar objects, masers, and other observational signposts of recent or ongoing star formation neai" 
the rims of GLIMPSE bubbles (N49 in Figure 2 provides an example, with two luminous young 
stellar objects and an ultracompact H II region visible on the lower rim of the bubble [17]). To 
date, most studies of triggered stai^ formation have focused on individual, "best-case" regions, very 
round bubbles with prominent sub-clusters or smaller bubbles on their rims [32-34]. In spite of 
the long-standing, widely popular idea that supemovae trigger star formation, far fewer candidate 
triggering sites have been identified near supernova remnants than near H II regions. 
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To establish that triggering merits investigation as an important mode of star formation as 
opposed to an astrophysical curiosity, the prevalence of triggering sites must be established sta- 
tistically, and in an unbiased fashion, among representative samples of H II regions. Both the 
GLIMPSE and MWP bubbles catalogs include flags for hierarchical structure, identifying smaller 
bubbles that could be the "daughters" of larger, "parent" bubbles. Among the MWP bubbles, 29% 
are members of hierarchies [27] . Triggered star formation need not produce observable daughter 
bubbles if the second generation of stars is too young to have produced H II regions. Thompson 
et al. [35] recently found a strong correlation between young stellar objects identified as luminous 
mid-IR point sources and the rims of bubbles from the Churchwell et al. [16] catalog. Corrobora- 
tion of this correlation using the more complete sample of bubbles from the MWP would greatly 
strengthen the case that triggering is a prevalent mode of Galactic star formation. 



5. Summary 

In this review, I have given a brief update on recent, multiwavelength observations, particu- 
larly in the mid-IR and X-rays, that have revealed wind-blown and radiation-dominated H II region 
structures. Observations corroborate the basic predictions of H II region theory but reveal impor- 
tant differences, too. Theory still struggles to explain both the existence of dust within energetic 
H II regions and its effects on H II region structure. X-rays from hot, wind-shocked plasma ai^e 
frequently observed to fill large, central cavities in giant H II regions, but the X-ray luminosity is 
more than an order of magnitude lower than predicted by wind-blown bubble theory. 

The size of the central cavities as a fraction of H II region volume appears to increase with 
increasing ionizing luminosity (Figures 2 and 3). This trend puts dust at larger distances from 
the ionizing stars in stai^burst regions, which reduces the effective temperature of the dust com- 
pared to static or thermal pressure-dominated H II region models. Single-band IR diagnostics of 
extragalactic stai^ formation rates (e.g. [36]) have become increasingly popular- in the era of large, 
high-resolution surveys from Spitzer and Herschel. Because the effective dust temperature sets the 
shape of the IR spectral energy distribution in unresolved regions, feedback must be taken into 
account when calibrating the IR diagnostics. 

Feedback in large, starburst regions like 30 Dor presents a particularly complicated problem. It 
would be premature to assume that a single source of feedback, radiation pressure (LI 1), dominates 
such regions until we achieve a better understanding of the contributions from massive star winds 
and supemovae. 

H II regions are among the most beautiful objects appearing in astronomical images. The 
strong aesthetic appeal of these images helped to motivate MWP volunteers, the vast majority of 
whom were not professional scientists, to spend tens of thousands of person-hours finding and 
measuring several thousand H II regions in the form of IR bubbles. The large MWP database [27] 
provides an unmatched resource for statistical studies of H II region evolution and star formation 
triggered by massive star feedback. I expect that the MWP will spawn many follow-up investiga- 
tions in the coming years, involving both professional researchers and citizen scientists. 
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